Geometrical aspects of a generalized statistical mechanics 
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We discuss here the use of generalized forms of entropy, taken as information measures, to char- 
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I. INTRODUCTION 

Our aim is to add a note on current efforts to inves- 
tigate the geometrical structure of the space of thermo- 
dynamic parameters for classical and quantum systems. 
This kind of studies has been developed by many authors 
(see, for instance, Refs. H S S S |I B S S @ ) , and has 
been applied to the ideal gas, the van der Waals gas, 
and magnetic systems that exhibit phase transitions. In 
such instances the Riemannian scalar curvature R of the 
parameters' space plays such an important role in de- 
tecting critical behavior that it has been interpreted as 
a measure of the stability of the thermodynamic system 
(typically, for non-interacting models one obtains R = 0, 
i.e. a flat geometry, while R diverges at the critical point 
for interacting systems). 

InRef.0 Wcinhold considered the thermodynamic sur- 
face given by the fundamental relation U = U({Xi}) in 
the (r + l)-dimensional Gibbs space (with coordinates la- 
belled by U,Xi, . . . , X r ), and obtained the components 
of the metric tensor of that space as the second deriva- 
tives of the internal energy U with respect to each pair 
of the r extensive parameters Xi . Ruppeiner Q focused 
attention on fluctuations of the thermodynamic magni- 
tudes and obtained the metric tensor via second moments 
of the fluctuations. The two ensuingrnetrics have been 
proven by Mrugala to be equivalent [3|. Another statis- 
tical path for reaching a Riemannian metric [ToL ITU IT^ 
in the space of thermodynamic parameters is that orig- 
inated in the works by Rao 0] and Amari (see, for in- 
stance, 01) m the field of statistical mathematics, and 
also by Ingarden 0, Janyszek |{|, and other authors in 
the field of thermodynamics and statistical mechanics. 
The concomitant information-theoretic approach is based 



on the concept of relative entropy and given in terms of 
the Boltzmann-Gibbs-Shannon entropy. 

There are also some recent efforts in the field of infor- 
mation geometry related with the generalized, nonexten- 
sive formulation of statistical mechanics [lol fig . Among 
these studies one finds the contributions by Amari, Na- 
gaoka and coworkers (0| among others) in connection 
with the geometrical structure in the manifold ofprob- 
ability distributions, and the analysis by Abe of 
the geometry of escort distributions, in connection with 
bit-variance and fluctuations of the crowding index of 
a multifractal. A new geometrical approach to thermo- 
statistical mechanics is introduced in 0], where the rel- 
evance of the approach within the contexts of nonexten- 
sive statistical thermodynamics is analyzed, showing that 
Riemannian geometry concepts yield a powerful tool. 
More recently, Naudts studies escort density oper- 
ators and generalized Fisher information measures. 

We will analyze here the geometrical approach to sta- 
tistical mechanics with regard to generalized measures of 
information, as the one usually denoted as Tsallis' nonex- 
tensive g-entropy 0, 
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for which we will show how to derive the metric tensor 
and other geometric quantities in the space of parameters 
that define the state of the system, p. Our considerations 
will refer to a generalized thermo-statistical context char- 
acterized by a given nonextensive index q and by a given 
prescription for the expectation values of relevant mag- 
nitudes. 

The paper is organized as follows: in Section [n] we 
summarize some basic results for the density operator in 
a generalized statistical framework. Section llTTl is devoted 
to exhibit some standard results on the geometrical as- 
pects of the space of parameters, while in Section Hvl ve 
present an extension of that study to a generalized statis- 
tical framework. For an arbitrary value of the nonexten- 
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sivity index q, we obtain the metric tensor of the parame- 
ters' space and discuss the meaning of various alternative 
ways of rewriting the metric tensor. Then, we calcu- 
late the scalar curvature in terms of the pseudo-partition 
function (given by Eq. (J2J below), the generalized mean 
values, and their derivatives with respect to the thermo- 
dynamic parameters. As an application we consider the 
ideal gas case in Section [V] Finally some conclusions are 
remarked in Section IVII 

II. THE DENSITY OPERATOR FOR 
GENERALIZED STATISTICS' OLM- VERSION 

The density operator that maximizes Tsallis nonexten- 
sive g-entropy Q for a given positive value of q is writ- 
ten, within the optimized Lagrange multipliers (OLM) 
formalism, as |2(j 

where mi , . . . , m r are the generalized expectation val- 
ues of r quantum operators. They are considered to be 
known as prior information and are given by the following 
generalized prescription [23, El 

Tr pip, 

mi = {Fi) q = g , i = 1, . . • , r (3) 

(for the sake of notational simplicity, we do not use along 
this article any label q for p nor m, despite the fact 
that they do depend on that index; we keep, however, 
the subindex in quantities such as S q and Z q ). The set 
of Lagrange multipliers {/3 1 , . . . , (3 r } is such that it fits 
the restrictions J2J in the procedure of constrained ex- 
tremization of S q when one is working within the OLM 
formalism po| . i.e. the restrictions are rephrased as cen- 
tered mean quantities, Tv[p q (Fi — m,)] = 0. This for- 
malism is an alternative to the TMP one of Ref. It 
has been established |2^,|2^| that, within OLM method, 
the correspond to the physical intensive parameters 
(e.g. the multiplier associated with the Hamiltonian op- 
erator, A = H, is {3 1 = (3 = l/(k B T) where T is the 
physical temperature of the system). In this framework, 
the pseudo-partition function 

Z 9 = Tre q (-^/?{Fi-mi)\ (4) 

is defined in such a form that the density operator is 
automatically normalized, i.e. Trp = 1. Notice that no 
Lagrange multiplier associated with this normalization 
condition for p is explicitly mentioned in writing down 
the equilibrium density matrix; instead, here we choose 
to introduce the pseudo-partition function in the expres- 
sion for p. In all these expressions, e q (x) stands for the 
q-exponential function 

e q (x) = [l + (l-q)x} 1 /^ (5) 



with [X]+ = 10(1), being O the Heaviside step func- 
tion. Finally, the "extensive limit" corresponds to the 
situation q — ► 1, and q — 1 becomes a measure of the 
degree of entropic nonextensivity. 

III. GEOMETRICAL ASPECTS OF 
PARAMETERS' SPACE IN STANDARD 
STATISTICS 

One of the key ideas for achieving our aims is to repre- 
sent the equilibrium density operator, and related quan- 
tities, in an r-dimcnsional space whose coordinate axes 
correspond to the parameters (i 1 , . . . , f3 r (called the sta- 
tistical temperatures). These are associated to the opera- 
tors Fx , . . . , F r , the mean values of which are assumed to 
constitute the known prior information. One can con- 
sider "neighbor" quantum states in the r-dimensional 
space of parameters and, from a notion of distance be- 
tween states, derive the metric tensor. Janyszek and 
other authors 0, El 01 IS HI have done so for the equi- 
librium density operator that maximizes von Ncumman- 
Shannon-Gibbs entropy, i.e. p = Z^ 1 exp(— (3 l Fi) 
with Z = Tr exp(— (3 l Fi). For the case of commut- 
ing operators F\, the metric tensor can be defined as 

9ij(0) = (djhxp dihxp) , i,j = l,...,r (6) 

where (...) = Tr (p . . .) and di = J—. For the case 
of non-commuting operators, this expression should at 
least be symmetrized. However a different derivation can 
be adopted, introducing the notion of information dis- 
tance between two quantum states, say p and cr, through 
the relative entropy or information gain (also known as 
Kullback-Leibler divergence), which is given by 

K(p\\a) = Tr pQnp - In. &) (7) 

The derivation of the metric tensor is as follows (see, 
for instance, 0): first a symmetrical (squared) infor- 
mation distance is constructed, as the sum D 2 (p,cr) = 
K(p\\a) +K(a\\p); this quantity is evaluated for a pair of 
equilibrium density operators, p({(3}) and p({/3 + A/3}), 
which arc neighbor states in the parameters' space. The 
ensuing expression is then expanded as a Taylor series 
for small displacements A/3 1 in all directions, resulting 
in vanishing zeroth and first order contributions. The 
non- vanishing second order terms thus allow for the defi- 
nition of the metric tensor from the local square distance 
dD 2 = gij d[3 l d(3 j as 

9iM = -(dA^P) (8) 
= djdilnZ (9) 

or, recalling that (Fi) = —di\nZ, 

9iM = -dj(F t ) (10) 

For classical systems or commuting operators, the def- 
initions for gij given by Eqs. (jSJ and (JSJ do coincide. In 
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these cases it is also seen that the fundamental tensor is 
related to thermodynamic fluctuations, as the component 
ij of the tensor corresponds to the covariance or second 
moment of the pair of operators Fi and Fj : 

9ij = ((Fi — — (Fj)) ) (11) 

Let us comment that for non-commuting operators an in- 
tegral expression for the covariances is introduced, based 
on this connection with the tensor g^ -see Ref. for 
details-. Finally, we notice that the metric tensor, be- 
ing constructed from second derivatives of the Kullback- 
Leibler divergence, corresponds to the Fisher information 
matrix. 

The metric tensor is one of the key ingredients in the 
geometric approach to thermostatistics. All the different 
viewpoints summarized here, Eqs. jB)|— l(TJ) , have been 
worked out by many authors in the framework of stan- 
dard (q — 1) statistics, and a number of applications to 
different physical models have been analyzed in this way. 
We present in the following section the formalism leading 
to analogous results for generalized statistics. 

IV. EXTENSION OF GEOMETRICAL 
CONCEPTS TO NONEXTENSIVE STATISTICAL 
FRAMEWORKS 

A. The fundamental tensor 

In the wake of the approach presented in Ref. |6| and 
working here along similar lines, we consider a general- 
ization of the definition of information gain to nonex- 
tensive set ting s of index q ^ f. The quantum q- 
divergence |24l I2M [2(| for two normalized operators p 
and a has been defined as 

KM\&) = T~~ (Trp-Tr/5^ 1 -") (12a) 
1 — q v ' 

generalizing the quantum Kullback-Leibler divergence or 
relative entropy, Eq. (JJJ. K q can also be written in the 
following equivalent form, 

K q (p\\a) = Tr^(ln 9 p - \n q a) (12b) 

1 — q_-i 

where the so-called (/-logarithm ln g x = - 1 _ for x > 
(excluding x = if q > 1), represents the inverse function 
of e q (x). It is straightforward to show that e q (x) and 
h\ q (x) tend to the exponential and logarithm functions, 
respectively, for q sufficiently near 1 (standard limit). 

The next step in our derivation is to extend the con- 
cept of (squared) information distance, D 2 (p,a). One 
can see that K q (p\\a) 7^ K q (a\\p), except for the par- 
ticular cases: (i) q — 1/2 (Hellinger distance, with the 
conditions Tr/5 = Tra = 1), and (ii) q — (since Kq = 
for any pair of normalized density operators). In general, 
then, for our purposes we make use of the symmetrized 
sum 

K q (p\\a) + K q (a\\p) (13) 



as information measure. From a Taylor expansion of this 
quantity around a certain set of values for the parameters 
{Pi} (i.e., a given point in the r-dimensional /3-space), 
that characterize the state of the system, we determine 
the form of the fundamental tensor in a generalized con- 
text (in fact, we determine a family of metrics which are 
characterized by the single index q, with q £ M + ). Let 
us then construct the first few terms of the associated 
Taylor series, taking p = p(j3) and a = p((3 + A/3): 

• The zeroth-order term is identically zero since 
K q (p\\a) — if and only if p = a. This was also 
the case for the standard (q = 1) context. 

• From Eqs. (|12l) and (|13l) . the first-order term van- 
ishes identically as in the standard case. This is 
a consequence of the shape adopted for the sym- 
metrized information gain. In other words, the ex- 
plicit form of the density operator does not play 
any role at this level. 

• From the second-order terms in our Taylor expan- 
sion of the symmetric information measure i|13fl . 
we obtain the following result, which we interpret 
as the generalized fundamental tensor within the 
nonextensive picture: 

9 (q) = 9Tr (p-^^l-) (14) 

= - qTi {pww) (15) 

where p{(3) is the OLM density matrix (J2J. We 
notice that this result can also be regarded as 
the ij element of the generalized Fisher informa- 
tion (GFI) matrix = (dij^j (see, for in- 
stance, an d references therein). Interestingly 
enough, from Eqs. ©, ©7 an d lfl5|l we gather 
that this quantity -when expressed as a trace- ex- 
hibits shape invariance, irrespective of the value of 
q. Indeed, the form of the GFI matrix is seen to 
be the same for any value of the nonextensivity 
index q, and it also coincides with that of the stan- 
dard case. Differences arise, of course, when one 
introduces explicitly the actual density operator, 
which may correspond either to the extensive or to 
the nonextensive scenario. 

Summing up, we have found the fundamental tensor of 
the space of parameters within a nonextensive statistical 
framework, in terms of first logarithmic derivatives of the 
generalized density operator, g$ = qTi (p dj In p di In p), 
as well as in terms of second derivatives of In p, Eq. 115fl , 
with p given as in Eq. J5J). Notice that, in both expres- 
sions, the standard logarithm function is involved. 

It is instructive to look at alternative formulations of 
the results already obtained. After a little algebra we 
can rewrite each of the expressions mentioned previously, 
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in terms of variations of the logarithm of the pseudo- 
partition function: 

<4 9) = q{-d jmi -(q- l)Z 2 ^Tv [p 2 ^ 1 X 

x(5 q Fj - dj lnZ q ){5 q F % - d t \nZ q )]} (16) 

and 

gjf = q{-8 3 In Z q di In Z q + Z^Tr [p 2 ^ 1 x 

x (S q Fj - dj In Z q ){S q Pi - di In Z q )]} (17) 

respectively, where the generalized deviation operators 
are defined as 5 q Fi = Ft — rm with rm given by Eq. J2J. 
Comparing these expressions we immediately obtain a 
useful (symmetric) relation for the first derivative of a 
given mean value with respect to any of the Lagrange 
parameters that writes 

djim = ditrij = dj In Z q d, In Z q - qZf^Tr [p 2 ^ 1 x 
x (SqFj - dj In Z q )(S q Fi - d l In Z q )] (18) 

from which the components of the metric tensor can fi- 
nally be cast in compact form as 

g\f = -a jmi -(q- l)d 3 In Z q d t In Z q (19) 

This covariant tensor of rank 2 is seen to be symmetric 
(by use of Eq. l(T%|) ). Its conjugate or reciprocal tensor, 
g( q > % i , is a symmetric tensor of contravariant valence two. 

Yet another form of representing the g-metric tensor 
involves the use of fluctuations, as in Eq. Ijlll) . This point 
is discussed with some detail in Ref.HjJ where generalized 
expressions for operator variance and covariance are also 
considered. 

We see from (|19fl that, in order to evaluate the gen- 
eralized metric tensor for a given system, we need the 
logarithmic derivatives of the pseudo-partition function 
with respect to the Lagrange parameters, as well as the 
behavior of the generalized expectation values of the rele- 
vant operators in /3-space. We end up by remarking that 



the standard limit for the fundamental tensor is correctly 
obtained in all the equivalent forms developed in this sec- 
tion. 



B. The scalar curvature 

Let us now address the computation of the curvature 
of the Riemannian space of parameters. We remind that 
in the extensive q = 1 framework Q the metric tensor 
G = (gij) corresponds to the matrix of second partial 
derivatives of the function /({/?}) = InZ with respect 
to every pair of parameters, i.e. the Hessian matrix of 
InZ, G = HilnZ) (see Eq. © ). As a consequence, the 
curvature tensor and the scalar curvature of the Riemann 
space, being constructed from derivatives of the elements 
of G, adopt a relatively simple form. Using Eq. ©, the 
results are the following [(|: 

Rijkl — -^g mn {f,milf,njk ~ f,mikf,njl) > f = hi Z 

(20) 

for the curvature tensor, where = dif = din Z/df3 l 
and so on, and 

R=-Ri2i2, g = \G\ (21) 
9 

for the scalar curvature, with g = \H(f)\. Equation l|21(l 
is valid in the particular case of considering only r = 2 
parameters. 

In a nonextensive setting characterized by a given fixed 
value of q, we arrive at more involved expressions that, 
however, reduce to the corresponding q — > 1 standard 
limit. For the sake of brevity, we define here 

f = \nZ q 

(using the standard logarithm) so that the Riemann- 
Christoffel tensor reads 0, ITlj 



Rifkt = -{Q - l){f,ikf,ji ~ f,af,jk) - 5 (<z);5t l-^(m s ,jim t ,ik - rn St jkm t ,u) + 

~ !) [fA m t,ikf,ji - m t ,uf,jk) + (f,ikm s ,ji - f,um s ,jk)f,t\ + (? - l) 2 {f,ikfji - f,uf,jk)f, s f,tj (22) 

I 



(with the notation m s ^j = djdim s ). This is one of our 
main results in the present effort. The m, correspond 
to generalized expectation values, and / stands for the 
logarithm of the pseudo-partition function. Apart from 
those g-dependent quantities, we appreciate in Eq. I|22|) 
the fact that (i) some new terms, of order (q— 1) and (g — 



l) 2 , emerge in the nonextensive framework with respect 
to the standard result l|20|l . and (ii) the metric tensor 
itself contains an explicit contribution with (q— 1), arising 
from Eq. I|19fl • It can be checked that the only terms that 
survive on letting q — > 1, exactly correspond to Eq. 1)20(1. 

We can now compute, by appropriate contractions, 
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first the Ricci tensor Rjf) = g( q ) ll R$ kl and then the 

scalar curvature = g^ q ^ k R^. A combinatorial 

analysis shows that, for an r-dimensional Riemann space, 
the number of independent components of the (4, 0)-type 
Riemann-Christoffel curvature tensor is r 2 (r 2 — 1)/12, 
due to symmetry reasons |l l| . At this point, in order to 
fix ideas and also for application purposes, we address the 
particular case of r — 2 Lagrange parameters character- 
izing the system. In this case it suffices merely to obtain 
an expression for R^iz and, then, the desired scalar cur- 
vature R^ q ' can be computed making use of a relation of 
the form given in Eq. I|21l) (notice that this equation is 
valid for any 2-dimensional metric). On the one hand, 
we calculate the 2 x 2-determinant g^ = |G^'|: 



3. the term 

+ f,t( m s,22f,ll 

times 



(g_l)l fl (9).t [f s ( mt 



g^^-_[f, s {m t ,nf,22 - 
2i/,i2)] equals (q 



m t ,i 2 /.2l) + 

- l)/(2 5 (9) ) 



ll 


h 


( 


In 


1,12 




/,21 /,22 


m 2 ,i 


m 2 , 2 




mi, 2i 


"T-1,22 




mi,u 7711,12 



ll 


h 


( 


In 


J, 12 




f,21 f,22 


mi,i 


n%i,2 




m 2 ,2i 


m 2 ,22 




m 2 ,U 7722,12 



26) 



,(«) - 



mi,i mi >2 
m 2 ,i m 2 , 2 

mi,i mi,2 
m 2 i m 22 



(q-r 

- (q - 1) x 

1,2 



f,i fa 

"22,1 7722,2 



/,l/,2 
2 
2 



/,! /,2 

mi,i mi 2 



(23) 



and, on the other hand, we compute separately each term 
appearing in Eq. for (ijkl) = (1212): 

1. the term -(q - l)(/,u/, 2 2 - /,i2/,2i) equals 



4. and the term -(<? - l) 2 (/,n/, 2 2 - /,i2/,2i) 
s * /, s /,i equals (g - l) 2 /.g (<z) times 



f,ll f,12 
f,21 f,22 



f,l f,2 
"22,1 ™2.2 



/,1 
mi i 



h 

mi 2 



(27) 



/,12 
/,21 /,22 



-(9-1) 

which is proportional to 



(24) 



2. the term — j g^ st (m Si2 2 m t,n — m s,2i^i,i2) equals 
l/{Ag^) times 



mi : i mi2 m 2 , 2 

7721,11 "11,12 "11,22 
7T2i,i2 772l,22 m-2,22 



+ (9-1) 



/I 



/,1 /,l/,2 J ,2 
7T2l,ll 7721,12 "11,22 
mi 12 772i 22 m 2 ,22 



(25) 



Thus ^1212 is obtained as the sum of the contributions 
in expressions l|24[l - (|27[) . with their corresponding pref- 
actors, and with g^> given by Eq (|23|l . Notice that we 
have arranged all the second order determinants as Ja- 
cobian ones, involving the two-parameters functions mi, 
7722, and / as well as their first derivatives with respect 
to both parameters. For brevity, we will denote the Ja- 
cobian determinant d((f>, -0)/9(/3 1 ,/3 2 ) simply as J(4>,ijj)- 
Also, the Hessian determinant of / is to be denoted as 
H(f) wherever it appears (see Eq. (|2"l|l ). Summing up, 
the generalized scalar curvature of 2D-space becomes 



( ff ( 9 ))2 



7721,1 772i,2 772 2 ,2 

7721.11 "21,12 ™1,22 

7721.12 7721,22 772 2 ,22 



+ (9-1) 



f,l /,l/,2 /,2 
7721,11 mi, 12 mi, 22 

mi. 12 mi. 22 m 2 ,22 



J(mi,m 2 )\H(f)\ + 



+ 2 [ J (f, m 2) (J(f.l,mi. 2 ) - J(/,2,772l,l)) - J(/,mi) (J(/,i,m 2 , 2 ) - J(f,2,m 2 ,i))] 



(28) 



with 

gM = 7(7721,7722) + (q - 1) [li J(/,m 2 ) - f, 2 J(/,mi)] 

(29) 

In these expressions it is not difficult to recognize the 
extensive limit, that coincides with the result given, for 



instance, in Ref.lff Finally, one could regroup all the con- 
tributions to R^ and analyze the ensuing departure from 
the extensive situation, which may eventually shed some 
light on the behavior of the system under study. Those 
terms that do not contain any explicit factor (q — 1)™, or 
"zeroth-order terms" , are clearly the ones that survive 
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in the limit q — > 1, while higher order terms in (q — 1) 
represent the successive corrections to the extensive case, 
when q is close to one. 



V. THE IDEAL GAS CASE 

We address here a classical ideal system consisting of N 
free particles in contact with a heat bath at temperature 
T and pressure P through a movable wall (the volume 
is variable). Considering the P-T, or Boguslavski, distri- 
bution, the standard partition function (usually denoted 
as Yn(P,T) in textbooks) is given by 



Z\(a,(3) 



1 



Mh 3N 



„-(JV+i) 



dV I d 3N qS N P e- aV ^ H 



/27rm\ 



3N/2 



(30) 



where H denotes the free particle Hamiltonian for N par- 
ticles of mass m, and the intensive parameters are taken 
to be a = P/kT and = 1/kT. 

One can employ Eq. (|3U[) as a starting point to cal- 
culate the nonextensive thermodynamical magnitudes of 
the system, by recourse to integral representations based 
on the definition of the Euler gamma function (see, for 
instance, the appendix in Ref. OMJ for details and a dis- 
cussion regarding its validity). Performing the so-called 
Hilhorst transform, which is applicable if the nonexten- 
sivity index q is greater than 1, the generalized pseudo- 
partition function can be written in the form 



Z q {a,p) = Zi{a,p) x 



(A 



\N- 1 



(A) 



fiV+l 



iN-l 



Oft) [i-(g-i)l,]^ 



(31) 



where for brevity we have defined the quantity I q = 
a(V) q + (3(H) q = aVq + PUq, which is a convenient com- 
bination of the intensive parameters and corresponding 
generalized mean values for volume and energy, respec- 
tively. The range of q for which Eq. I|31|) is valid is given 
by < g — 1 < (|iV + l) _1 . Performing the same sort of 
real integral representation as before, one can compute 
the ^-expectation values of the relevant operators. The 
following interesting result ensues 



V a 



N + 1 



and 



U a 



3N 



(32) 



while I q = |iV + 1. It is remarkable that these results 
are independent of the value of the index q (as long as it 
belongs to the range of validity of the integral representa- 
tion employed). Also, for later purposes, we notice that 
U q does not depend on the pressure, while V q is a function 
of a only. In both expressions, N is the (fixed) number 
of particles that constitute the system under study; note 
that in the "equation of state" there is a factor N + 1 



instead of N as we have not taken the thermodynamic 
limit which allows for that replacement. 

With these ingredients we can compute the generalized 
metric tensor and other geometrical quantities of interest 
defined in Sectional From Eqs. J3U and l|3*2*)l. for the q- 
metric tensor Ijl9|l we obtain the following 2x2 symmetric 
matrix 



4? 



JV+l 



_ ot 





37V 
2iP 



"(9-1) 



m 

3N(N+1) 
2a/3 



for 1< q < 1 + (fiV 



3N(N+1) 
2a/3 

(«0 2 

(33) 

1 + gives the 



l)" 1 . The limit q 
correct diagonal result already obtained by other authors 
(Refs.0and|3, for instance). Here the changes introduced 
by the nonextensive treatment appear as contributions 
of order {q — 1): a dependence with N 2 appears since 
we have not considered thermodynamic magnitudes per 
particle (as Brody and Rivier, for instance, do). The 
2 x 2-dctcrminant of the matrix l|33l) acquires a simple 
form in terms of the intensive parameters a and /3, and 
the number of particles N: 

We compute explicitly the expressions in l)24 |l — 1|27| ) (with 
their corresponding prefactors) and arrive, respectively, 
at 

2 



(g-1) (3N(N + 1) 



{-l + (q-l)(-N+l)} 




(g-1) ( 3N(N+1) 

(g-1) ( 3N(N+iy 2 
gii) \ 2a 2 /3 2 



X {-( g -l)(-AT+l)} (35) 



Adding these four terms and multiplying by the factor 
we are led to R^ q \ The computation of the gen- 
eralized scalar curvature in this example yields the out- 
come: R( q ) = 0. It is noticeable that this result (i) is 
independent of the intensive parameters a and (3, and 
(ii) is also independent of the nonextensivity index q. 
We can then assert that in this case the outcome for 
the curvature in our geometrical approach, performed 
within the framework of a nonextensive statistical set- 
ting, agrees with the developments obtained by different 
authors within the structures of the extensive formalism 
-see, for instance, Refs. 0, Hi-. 



VI. CONCLUDING REMARKS 

We have discussed geometrical aspects of the space 
of parameters for a physical system, within the frame- 
work of a nonextensive statistical setting of fixed index 



7 



q. For that purpose, we have employed a generalization of 
the quantum Kullback-Leibler divergence as information 
gain. In particular, we have derived the pertinent gen- 
eralized metric tensor which is connected with Fisher's 
information matrix. Our main findings in the present ef- 
fort are represented by Eqs. (|19J) and Q28|l . which pave the 
way for deriving the metric and curvature tensors from 
a generalized partition function for a fixed value of the 
nonextensivity index. 

The studies addressed here are relevant in the investi- 
gation of critical behavior and phase transitions in ther- 
modynamic systems. Work on such sort of applications 
for interacting systems, using generalized information 
measures, is currently in progress. 

Finally, we remark that the use of alternative mea- 
sures of distance between pairs of density operators in the 
space of thermodynamic parameters could yield an alter- 
native method -perhaps a more sensible one- to detect 



thermodynamic changes in certain physical systems ap- 
plying the geometric machinery discussed in the present 
effort. To mention just one such possibility, the Jensen- 
Shannon (JS) divergence can be used to give a measure of 
distance between two probability distributions. A gener- 
alization of the JS divergence in the framework of Tsallis 
statistical mechanics has been investigated by Lamberti 
and co-workers 31, 32], who obtained a monoparametric 
family of metrics. This and other alternatives will be the 
subject of future work. 



Acknowledgments 

The authors acknowledge financial support from CO- 
NICET and UNLP, Argentina. MP also acknowledges 
ANPCyT (PICT No. 03-11903/2002), Argentina. 



[2] 

[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9] 

[10] 
[11] 

[12] 



[13] 
[14] 



[15] 
[16] 



F. Weinhold, J. Chem. Phys. 63, 2479, 2484, 2488 (1975); 
65, 559 (1976) [17 

G. Ruppeiner, Phys. Rev. A 20, 1608 (1979); Rev. Mod. [18 
Phys. 67, 605 (1995); Erratum: 68, 313 (1996) [19 
R. Mrugala, Physica 125A, 631 (1984) [20 
R.S. Ingarden, Tensor N.S. 30, 201 (1976) 

H. Janyszek, Rep. Math. Phys. 24, 1, 11 (1986) [21 
H. Janyszek and R. Mrugala, Phys. Rev. A 39, 6515 
(1989) [22 
H. Janyzsek, J. Phys. A 23, 477 (1990) 
D. Brody and N. Rivier, Phys. Rev. E 51, 1006 (1995) [23 
W. Janke, D.A. Johnston and R. Kenna, Phys. Rev. E [24 
67, 046106 (2003); Physica A 336, 181 (2004) 

P. Le Corbeiller, Sci. Am. 191, No. 5, 80 (1954); [25 
L.A. Santalo, Vectores y Tensores con sus aplicaciones [26 
(EUDEBA, Buenos Aires, 1981), 12a. ed. [27 
David Poole, Algebra lineal: Una introduction moderna 
(Linear Algebra: A Modern Introduction) (Thomson, [28 
Mexico, 2004) [29 
C.R. Rao, Bull. Calcutta Math. Soc. 37, 81 (1945) 
S. Amari, Differential-Geometrical Methods in Statis- [30 
tics, Lecture Notes in Statistics 28 (Springer, New York, 
1985); S. Amari and H. Nagaoka, Methods of Information [31 
Geometry (AMS and Oxford University Press, 2000) 
C. Tsallis, J. Stat. Phys. 52, 479 (1988) [32 
A periodically updated bibliography on nonexten- 
sive thermostatistics can be found in the URL 



http: //tsallis . cat . cbpf .br/biblio .htm 

S. Abe, Phys. Rev. E 68, 031101 (2003) 

R. Trasarti-Battistoni, cond- mat/0203536 

J. Naudts, Open Sys. & Information Dyn. 12, 13 (2005) 

S. Martinez, F. Nicolas, F. Pennini and A. Plastino, 

Physica A 286, 489 (2000) 

C. Tsallis, R.S. Mendes and A.R. Plastino, Physica A 
261, 534 (1998) 

S. Abe, S. Martinez, F. Pennini and A. Plastino, Phys. 

Lett. A 281, 126 (2001) 

S. Abe, Physica A 300, 417 (2001) 

P.W. Lamberti, M.T. Martin, A. Plastino and O.A. 

Rosso, Physica A 334, 119 (2004) 

S. Abe, Phys. Rev. A 68, 032302 (2003) 

S. Abe, Physica A 344, 359 (2004) 

M. Tegmark, A. Taylor and A. Heavens, Astrophys. J. 
480, 22 (1997) 

S. Martinez, private communication (2004) 

M. Portesi, A. Plastino and F. Pennini, 

cond-mat/0510434 

S. Martinez, F. Pennini, A. Plastino and M. Portesi, 
Physica A 332, 230 (2004) 

P.W. Lamberti and A.P. Majtey, Physica A 329, 81 
(2003) 

A.P. Majtey, P.W. Lamberti and A. Plastino, Physica A 
344, 547 (2004) 



